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Abstract 

In a previous paper, we have considered the weakly nonlinear interaction of a pair 
of axisymmetric lower branch Tollmien-Schlichting instabilities in cylindrical supersonic 
flows. Here the possibility that nonaxisymmetric modes might also exist is investigated. 
In fact it is found that such modes do exist and, on the basis of linear theory, it appears 
that these modes are the most important. The nonaxisymmetric modes are found to exist 
for flows around cylinders with nondimensional radius a less than some critical value a e . 
This critical value a e is found to increase monotonically with the azimuthal wavenumber n 
of the disturbance and it is found that unstable modes always occur in pairs. We show that 
in general, instability in the form of lower branch Tollmien-Schlichting waves will occur 
first for nonaxisymmetric modes and that in the unstable regime, the largest growth rates 
correspond to the latter modes. 


This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NASl-18107 while the authors were in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 


i 





1. Introduction 

Recent interest in the development of supersonic and hypersonic vehicles has stim- 
ulated research into the different instability mechanisms which might cause transition 
to turbulence in high speed compressible flows. Apart from Tollmien-Schlichting waves, 
Gortler, and crossflow vortices, the possibility of inviscid instability modes must be taken 
into account for compressible flows. Here an investigation of the role of nonaxisymmetric 
lower branch Tollmien-Schlichting waves in supersonic flows with an axis of symmetry is 
investigated. 

Following the work of Smith (1979a, b) it is now well known that triple-deck theory 
provides a self-consistent asymptotic framework for the description of lower branch linear 
and nonlinear, two and three dimensional Tollmien-Schlichting waves. Thus it is known 
that finite amplitude Tollmien-Schlichting waves are stabilized by nonlinear effects as they 
cross the lower branch of the neutral curve. At higher disturbance amplitudes further 
downstream Smith and Burggraf (1985) have shown that a hierarchy of fully nonlinear 
states can be achieved. The planar compressible problem has been investigated using 
triple-deck theory by Smith (1988). At subsonic speeds lower branch disturbances are 
essentially unchanged from their incompressible forms. However, at supersonic speeds 
only three-dimensional modes can be unstable and the critical angle of propagation above 
which instability can occur increases with the Mach number. At hypersonic speeds Smith 
(1988) has shown that the modes are then fully nonparallel and a quasi-parallel theory 
based on triple deck theory fails. It can be shown that in this regime the modes have a 
structure similar to that which describes Gortler vortices in incompressible flows (see for 
example Hall (1983)). 

In a previous paper, Duck and Hall (1988), we investigated the linear and weakly non- 
linear theory of lower branch axisymmetric Tollmien-Schlichting instabilities in supersonic 
cylindrical flows. We found that such modes exist in pairs and that at a given Mach num- 
ber they occur only for a body radius less than a critical value. In the limit of small body 
radius both modes have wavelength going to zero with respect to the usual triple deck 
streamwise length scale. In the weakly nonlinear stage it was shown that, dependent on 
the input frequency, either mode can lead to a stable finite amplitude state. 

Here we generalize the above calculation to see if nonaxisymmetric modes can be impor- 
tant in the linear regime. We perturb an axisymmetric supersonic flow to three-dimensional 
disturbances with wavenumbers a and n in the streamwise and azimuthal directions. We 
show that at all supersonic speeds there is a finite band of non-dimensional body radii 
which can support two three-dimensional modes of the same azimuthal wavenumber and 
that in the limit of bod radius tending to zero one mode has wavelength tending to zero 
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whilst the other tends to infinity. Thus for these cylinders Tollmien-Schlichting waves with 
both short and long wavelengths compared to usual streamwise triple deck scale can occur. 
At higher azimuthal wavenumbers and larger body radii we show how the modes relate to 
the planar three dimensional modes of Smith (1988). We shall show that, on the basis of 
linear theory, it is the three-dimensional modes which are the most dangerous since they 
occur at the lowest Reynolds number and when they occur have the largest growth rates. 

The procedure adopted in the rest of this paper is as follows: in §2 we formulate the 
linear stability problem for supersonic cylindrical flows. In §3 we determine the eigenrela- 
tions for three-dimensional TS waves around a cylinder with radius comparable with the 
upper deck thickness. In §4 the eigenrelation is investigated in the limit of small and large 
cylinders and also in the high azimuthal wavenumber limit. Finally in §5 we describe some 
results and draw some conclusions. 

2. Formulation 

We shall be concerned with the linear stability of an axisymmetric boundary layer on a 
cylindrical body of radius a*, in a uniform supersonic stream of velocity U ^ aligned with 
the axis of the cylinder. 

If L* denotes a typical streamwise length scale (for example the distance from some 
leading edge), the kinematic viscosity of the fluid in the far field, then the Reynolds 
number Re is defined to be 

Re = Tr oo V/v' oa . (2.1) 

It is found useful to introduce a related parameter £, 

e = Re~*. (2.2) 

In this paper the Reynolds number is taken to be large, implying e is small. 

It is also assumed that 

S=^ = °( 1 ). (2-3) 

denotes the scale of the radius of the body. This follows the scale chosen by Duck and 
Hall (1988), Kluwick et al. (1984), and one of the scales chosen by Duck (1984), amongst 
others; this turns out to be an important choice of body radius scale, with curvature terms 
playing a crucial role in the physics of the problem. 

The problem then takes on a triple deck structure. The following non-dimensional 
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variables are defined 


X = 


z*-£* 

« 3 L* 


r = 


r 

e 3 L* ’ 


* — P-P 


V = 


pSo^, 


u 

(/• 


^ J M — Jf* > 


(2.4) 


ty = 


m 

t/* > 


c = 


e 


ZZJL 

t 2 L’ ' 


t = 

Here x*,r* and 6 are taken to be the streamwise, radial and azimuthal co-ordinates respec- 
tively at some suitable reference point on the body, and u*,v*,w* are the corresponding 
velocity components, c* denotes the speed of sound, p ^ the fluid density far from the body, 
p* the pressure, p ^ the far field pressure, and t* the time. 

We confine our study to the stability of the flow at a location on the body where the 
boundary layer thickness is 0(e 4 L*), which is thin compared with the radius of the body. 
The skin friction of the undisturbed boundary layer is then taken to be U*Xe 4 /L*, where 
A is some order one parameter. The flow is taken to be parallel, a completely rational 
approximation in this context because of the small, 0 (e 3 Z»*), streamwise length scale under 
consideration (although non-parallel effects could become important in a non-linear study). 

We consider first the upper deck of the triple deck, where r = 0(1). The solution 
develops as follows 

= 1 + e 2 u 1 (X,r,0,t)+, 


u 


v = + 


w = e 2 toi (X,r,0,t) + 


p = e 2 pi(X,r,6,t) + 


(2.5) 


c = + e 2 c 1 (X,r,0,t) -\ . 

is the Mach number of the external flow. The flow in this layer turns out to be 
completely irrotational, and may be reduced to the Prandtl Glauert equation (expressed 
in cylindrical polar co-ordinates), for the pressure pi. 


(1 — M^pm + -pir + pir r + —pieo = 0. 


( 2 . 6 ) 


The solution of this will be deferred, until the appropriate boundary conditions have been 
ascertained. We consider next the main deck, which corresponds to the transverse scale 


y = T ~FU~ = °* 1 )' 


(2.7) 
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Here curvature no longer explicitly plays a leading role, and the solution takes the form 

a = U o (y) + eA(X,0,t)U! i (y) + -" t 


v = -e 2 Ax{X,0,t)U o (y) + ---, 

w = e 2 D(X i 0,t)/Uo(y) + • • • , (2.8) 

P — Ro[y) + cA(X,0, t)Ro u (y) + • • • , 
p = e 2 P(X,0,t) + 

where U 0 (y) and Ro(y) are the axial velocity and density distributions (radially), corre- 
sponding to the undisturbed boundary layer. The function D can be written explicitly 
as 

D x = -?L. (2.9) 

The limit as y — * oo of the transverse velocity component yields an inner boundary condi- 
tion for the upper deck, viz 

Pir |r=a = AxX’ (2.10) 

As y — ► 0, (2.8) violates the no-slip condition, and this is subsequently rectified by the 
inclusion of the third layer, the lower deck, where 

Y = yje = 0(1), (2.11) 

and 

u = eU{X,Y,0,t), 
v = e*V{X,Y,0,t), 

w = eW(X, Y, 0,t), (2.12) 

p = e 2 P(X,0,t ), 

P — Ro{0) + epi(X,Y ,0,t). 

It is now possible to scale out a number of the physical constants. Following Kluwick 


! 

I 
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et al. (1984) and Duck and Hall (1988), this is achieved as follows 

x = 

y = c’lA-l(r U) /r 00 )ly, 

p = CiA?P, 

U = C*\*(T w /Too)»U, 

V = cIAt^/Too)^, (2.13) 

W = c^x^t^t^w, 

A = cIa-i^/t^Ia, 

a = dX-iiTv/T^a, 

t = CiX-ziTv/T^t. 

Here C is the Chapman constant arising from the linear viscosity law 


/P//C = C{T/T W ), 


(2.14) 


T w being the non-dimensional wall temperature. The governing equations in the lower deck 


then turn out to be a form of the three-dimensional unsteady boundary layer equations, 


namely 

dU , j T dV , V du , WdU _ d 7 U _ dP 

at ^ u ax t v dY ' r a ae ~ W 7 ax' 


aw , rjdw , v aw , w aw _ a 3 w _ i ap 
at ' u ax ^ v ay ' a ae ~ aY * a ae > 

au I av , i aw _ n 
ax by "*■ a ae ~ u > 

with boundary conditions 


(2.15) 


u = V = W -0 on y = 0, 


U —* Y + A(X,0,t) as y — *■ oo, (2.16) 

W -> D/Y as y oo, 

where 

Dx = — Pe • (2-17) 

a 
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If we write the solution of (2.6) subject to (2.10) symbolically as 

P = Z{A} (2.18) 

then the problem is effectively closed. The system described above is inherently non-linear, 
and as such generally requires a numerical solution; however we now seek solutions which 
are small perturbations about the undisturbed state, namely 

U = Y, V = W = P = A = 0; (2.19) 

this amounts to a linearized stability analysis of the boundary layer, and forms the basis 
for the following section. 

3. Linear stability 

Here the undisturbed flow is perturbed by a small amount h « 1, and we seek to 
investigate the growth or decay of these disturbances. We seek a solution of the form 

U = Y + hUi -I , 

V = hV i + 

W = Wi + -, (3.1) 

P = hPi , 

A = hA\ + • • • , 

in particular the leading order perturbation terms are assumed to take on the following 
particular form 

lh = ui(Y)Ei + c.c., 

V I = t>i {Y)E l + c.e. t 

Wi = tbi[Y)Ei + c.c., (3.2) 

Pi = P\Ei + c.c., 

A 

A\ = AiE\ + c.c., 

where 

Ei = exp{z[aX + nQ — Ot]}, (3.3) 

and “c.c.” denotes a complex conjugate. Pi,Ai,a,£l are in general complex constants; n 
is the (integral) azimuthal wavenumber. 
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The problem at 0(h) reduces to the following linear system 

Ui{— if! + iaY} + t)i = taPi + tiiyy, 

thl'f — Oh + iaY} = ~~P\ + Uq yy, (3.4) 

iau x + viy + = 0. 

The solution to this system is best achieved by introducing the variable 



a a n „ 

a = Uiy H ttqy. 

(3.5) 


ota 

Utilizing (3.4), the following simplified system results 



oyy + [*n — iaY ] a = 0; 

(3.6) 

the solution (which remains bounded as Y — ► oo) is 



a = 

(3.7) 

where 

. (to:) ^fl .. . i_, 

i = — + M*y. 

(3.8) 

Invoking, further, the boundary conditions as Y — ► oo demands 



l f°° * 

(ia) *B / Ai(()d( = A u 

J £o 

(3.9) 

where 

(ta)?n 

Co 

a 

(3.10) 

Evaluating the first two equations of (3.4) on Y = 0, and combining yields 



(*a)ajBilt'(Co) = (tad- 

(3.11) 

To complete this solution we require the solution from the upper deck. We have 


' d 2 Id n 2 2 ft , , 2 \1- o 

[dr 2 + r dr r 2 “ ^ M °° )J Pl “ °» 

(3.12) 

where 

Pi(r) = C7Ai[px(r)Ei + c.c.], 

r = clAl(r./roo)lr, 

(3.13) 

subject to 

Plr | r=a= -Ot*A U 

(3.14) 
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and 

pi — > 0 as r — > oo. 

The solution to (3.12), subject to (3.14), (3.15) is 


- _ iaAiKnliary/M*, - 1 ] 

Pl r ~ y/Ml-lKUiaay/Ml -1] 


(3.15) 


(3.16) 


(assuming > 1). Finally, a non-trivial solution to (3.9), (3.11), (3.16) exists only if 
the following (dispersion) relation is satisfied 

4(6>) _ M*! 1 + &]Knliaay/Ml ~ *1 

M*W y/Ml - 1 K' n [iaay/Mi - 1] 

Notice that setting n = 0 (axisymmetric mode), and defining 

a = (M^-l)la, 


(3.17) 


H = (A£-l)ta, 
a = (M^ - l)“ia, 


(3.18) 


6 — 


is)in 


yields the following dispersion relationship 

At'do) _ (ta)gjiro[taa] 
Ki[iaa] 


(3.19) 


is Mm 

which is the same dispersion relationship found by Duck and Hall (1988) (unlike the non- 
axisymmetric situation, the Mach number may be scaled from the axisymmetric problem). 

In the following section, we describe various asymptotic solutions of the dispersion 
relationship and discuss the numerical solution of (3.17) in §5. 


4. The dispersion relationship for a « 1 and a >> 1 

Here we shall discuss the limiting forms of the dispersion relationship in the limits of 
either very thin or thick cylinders. We recall that axisymmetric modes always occur in 
pairs and exist for a less than some critical value, with branches having a, f 1 tending to 
infinity when a — * 0. Suppose then that a — ► 0 in (3.17) in such a way that oca — *■ 0. From 
the series expansions of K n , K' n we deduce that (3.17) then reduces to 


Ai'(Zo) 

X 


r- \i n . 

M 3 - + 

a 
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where 



so that for neutral modes we must have that 

£o ~ 2.298^ , Ai'(Zo)/x ~ 1.001* s. 


and the neutral values of a, fi then become 


a 

0 




2 >2 98{h22l2.} 2 + 


(4.1a, b ) 


Thus unlike the two dimensional problem there are neutral solutions for a « 1 which 
have a, VI — > 0. This is a result of some significance because it means that for thin cylinders 
non-axisymmetric modes are the most important since they will occur at lower values of 
the Reynolds number than do the axisymmetric modes. 

In fact there is another asymptotic solution of (3.17) available in the limit a » 1. 
The second mode again has aa « 1 but | £o | now tends to infinity. For large values of 
| £o | we can replace dills! by — £ 0 {l — £ 0 2 + ■ • •} and we again approximate using 

the series expansion of the modified Bessel function. After equating the dominant real and 
imaginary parts of (3.17) we then deduce that 


a 


{ 


2’ >+ £n!n--l! \ 2 "+7 


a 


4r>— 3 ^ 

4n+l / . • • 


n 



(4.2a, 6) 


Thus a ~ a - '*'>+ 1 ^ in this limit and since increases monotonically from | to 1 when 
n — > oo it follows that the largest values of a correspond to n tending to infinity. It follows 
that the distances between the upper and lower branches of (3.17) for a — *• 0 increases with 
n. Thus the band of unstable wavenumbers increases with the azimuthal wavenumber n; 
this does not of course necessarily mean that the maximum growth rate will occur for 
n » 1. We will return to the latter point in the next section. 

Now we consider the structure of (3.17) in the limit a — * oo. Firstly we note that if we 
let a — * oo with n fixed then, if | aa |>> 1, (3.17) reduces to 

^■'(&>) (*°^) 

X («£-!)* 


which has no neutral solutions. Thus we must instead consider the double limit n — > 
oo,a — > oo but with ^ held fixed. This means that the wavelength in the azimuthal 
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direction is comparable with that in the streamwise direction. The asymptotic form for a 
modified Bessel function of large argument and order yields 


K„(iaa^/(Ml l - 1) _ -t'qay'Mj, - 1 ^ 

- 1) " b{ 1 - ^[Mi - Tp + ” ' ’ 


where we have assumed that M — l] < 1. It follows that neutral solutions of (3.1) 
again occur when £o — —2.298*3' so that 


1 . 001(1 - 


a 2 a 2 


n L 


7 

r , ,,.1 OL^a - n 

[ML- l]} 3 = — -{1 + 


n 


a? a 2 


} + 


or if we write /? = n/a we obtain 


1-001(1 - - 1]}* = aJ/r^i + + • • • , 


(5.3) 


which is the supersonic planar neutral dispersion relationship of Smith (1988). 

It is an easy matter to show that (3.17) has no neutral solutions for n,a » 1 with 
M — 1] > 1. Thus neutral disturbances exist only for azimuthal wavenumbers satis- 
fying 

n 2 > a 2 a 2 - l}. (5.4) 

We shall discuss the above asymptotic results after describing our numerical results for 
the dispersion relationship (3.17). 


5. Results and Discussion 

The dispersion relationship (3.17) was solved for the neutral state ( Imag{a } = Imag{U} 
0), for prescribed Moo,n, using a straightforward Newton iterative scheme. Results for 
Mqo = \/2 are shown in Figs la (0 — a curve), lb (a — a curve), for M = 5 in Figs 2a 
(fl — a curve), 2b (a — a curve), for Moo = 10 in Figs 3a (fi — a curve), 3b (a — a curve). 

As with the axisymmetric case considered by Duck and Hall (1988), for fixed a less 
than a c , the critical body radius above which the flow is stable) two possible modes exist. 
However, there is an important difference between this and the axisymmetric case, namely 
the behavior of the “lower branch,” which here has 0 — > 0, a — > 0 as a — > 0, whilst in the 
axisymmetric case it is found that U — >oo,a— > oo as a - >0. 

Further observations regarding these n ^ 0 results are 

1. as n increases, a c increases (indeed, all the a c determined here were greater than the 
a c obtained in the n = 0 case) . 

2. as Moo increases, for n fixed, a c decreases ; 
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3. kinks were consistently observed on the upper branch of the a — a curve for the 
smaller value of n. 

In Figures 1, we have also shown (as broken curves) the lower and upper branches 
predicted by the asymptotic analysis of §4. We see that the asymptotic theory accurately 
predicts the neutral wavenumbers and frequencies over a wide range of values of a. A 
similar agreement is found for the higher Mach number cases but is not shown in the 
figures. 

In Fig. 4, we have shown the result of some non-neutral calculations corresponding to 
Moo = \/2 and a = 2. The frequency fl is taken to be real and when the corresponding 
complex value of the wavenumber a is calculated we see that the curves are unstable 
for frequencies between the two neutral values where the latter exist. Furthermore, the 
growth rates initially increase with n so that, at least at the values of a and chosen, 
the non-axisymmetric modes become progressively more unstable with increasing n. In 
fact the most dangerous mode, i.e., the one with the largest growth rate occurs for n = 6 
for a = 2, Mqo = \/2. A similar trend was found at higher Mach numbers and in general 
we found that the value of n for the most dangerous mode increased with M*. 
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Fig. 2b. 
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